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Abstract. Improved limits as well as tentative claims for dark matter annihilation into 
7-ray lines have been presented recently. We study the direct detection cross section induced 
from dark matter annihilation into two photons in a model-independent fashion, assuming no 
additional couplings between dark matter and nuclei. We find a striking non-standard recoil 
spectrum due to different destructively interfering contributions to the dark matter nucleus 
scattering cross section. While in the case of s-wave annihilation the current sensitivity of 
direct detection experiments is insufficient to compete with indirect detection searches, for p- 
wave annihilation the constraints from direct searches are comparable. This will allow to test 
dark matter scenarios with p-wave annihilation that predict a large di-photon annihilation 
cross section in the next generation of experiments. 
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1 Introduction 



Gravitational effects on astrophysical scales give convincing evidence for the presence of dark 
matter (DM) in our Universe, an observation which is strongly supported by measurements of 
the cosmic microwave background anisotropies [1]. While the existence of DM seems firmly 
established, very little is known about the properties of the DM particle(s). In order to 
identify their nature, several complimentary search strategies are currently being employed, 
which fall into three classes: direct detection in shielded underground detectors, indirect 
detection with satellites, balloons and ground based telescopes looking for DM annihilation 
signals and DM production at colliders such as the LHC. 

In general, indirect detection experiments looking for DM annihilations have to deal 
with astrophysical backgrounds, which make it difficult to unambiguously extract DM signals. 
One exception is annihilation into mono-energetic photons, which would provide a striking 
evidence for DM in our Galaxy. Consequently, 7-ray lines are among the most important 
search channels for the indirect detection of DM [2]. Recently, there have been hints at 
a tentative 7-ray line at E-y ~ 130 GeV associated with a rather large annihilation cross 
section of (a x v ve \) ~ 1.3 • 10~ 27 cm 3 s _1 [3, 4] (assuming an Einasto profile, see subsequently 
also [5, 6]). Such a 7-ray line has been predicted by various authors (see e.g. [7-10]). This 
finding, although in slight tension with the upper limit (<J x v re \} < 1.0-10 - cm 3 s _1 set by the 
Fermi LAT Collaboration [2], has triggered a noticeable amount of theoretical studies [11-26] 
most of which aim at explaining the possible 7-ray signal in terms of new physics. 

In the following we show that in view of the impressive sensitivity of current (and 
upcoming) direct detection experiments, relevant constraints on the nature of DM can arise 
from these searches even if the interactions of DM with quarks and gluons are loop suppressed. 
In particular, we will explore the constraints on the annihilation cross section of DM into 
photons that arise from direct detection bounds. In order to keep our discussion as general as 
possible, we will work in an effective field theory (EFT) framework obtained by integrating out 
heavy degrees of freedom above a certain high-energy scale. Our model-independent analysis 
includes the most relevant effective operators of dimension up to 7, containing bilinears 
constructed from scalar and fermion DM fields which lead to annihilations into 77. 

We find that in the case of annihilation via s-wave, the constraining power of present 
direct detection searches is far below that of indirect detection experiments. For p-wave 
annihilation, on the other hand, the constraints from direct and indirect searches are com- 
petitive. Given the expected improvement of direct detection experiments, DM models with 
p-wave annihilation and large di-photon cross section will hence be testable in the near future. 
This finding underscores the possible complementary between direct and indirect searches in 
unraveling the precise nature of DM. 

Our work is organised as follows. After listing the relevant effective operators that can 
give rise to a mono-energetic 7-ray line in Section 2, we calculate in Section 3 the correspond- 
ing annihilation cross sections into a pair of photons. Section 4 is devoted to a comprehensive 
discussion of the loop-induced effective interactions relevant for DM direct detection. In this 
section we also estimate the associated nuclear matrix elements, assuming the dominance of 
a single operator. Applying recent bounds from DM direct detection experiments, we finally 
obtain in Section 5 lower limits on the suppression scale of the effective operators, which in 
turn translate into model-independent upper bounds on the 77 annihilation cross section. A 
summary of our main results is presented in Section 6. A series of appendices contains useful 
details concerning technical aspects of our calculations. 
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2 Operator analysis 



Throughout our analysis we will assume that both DM annihilation and direct searches can 
be described in terms of effective operators that are generated by integrating out heavy 
degrees of freedom. While such an EFT provides an excellent description of the low-energy 
processes involved in direct detection, the mass scale M* of new physics should be sufficiently 
high for the framework to be applicable for DM annihilation. In particular, scenarios with 
resonant s-channel annihilation are not covered by our analysis. 

In our work we consider only effective operators that lead to the process XX ~^ 77 > 
where x denotes the DM particle. As pointed out in [18], dimension-4 operators correspond 
to milli-charged DM, which is strongly constrained experimentally and will hence not be 
discussed here. If DM is a Dirac fermion D, the leading contributions then arise from 
dimension-5 operators of electric or magnetic dipole type. Such operators give rise to long- 
range interactions between DM and nucleons at tree level and therefore to a sizeable direct 
detection cross section. The resulting constraints have been studied in [27]. The authors 
point out that DM with dipole-type couplings loses energy as it penetrates the Earth, so the 
strongest bounds arise from direct detection experiments in shallow sites. Nevertheless, these 
bounds still constrain the annihilation cross section into 77 to many orders of magnitude 
below the thermal cross section. An observable 7-ray line can hence not arise from such 
operators and we will therefore not consider these interactions any further. The same line 
of reasoning applies to other higher-dimensional operators that induce a tree-level coupling 
of DM to nucleons like the charge radius operator for complex scalar DM. Consequently, we 
will restrict ourselves to operators that are bilinear in both the DM and the photon fields. 

For real scalar DM R, we obtain two different dimension-6 operators 

q r = C R R 2 F^F^ U , Of = C?R 2 F^F^ , (2.1) 

where F^ v = d fM A u — d u A fl denotes the usual electromagnetic field strength tensor, while 
pfiu _ i e fj.fpup^ with e 0123 = 1 is its dual. The corresponding operators for complex scalar 
DM C, called O and Op, are obtained from the above by replacing R 2 with C^C. Notice 

R C 

that due to dimensional reasons the Wilson coefficients C,l appearing in (2.1) all scale as 

For fermionic DM, the relevant operators are of dimension 7. In the case of Majorana 
fields M, we obtain two different sets of operators depending on whether the interactions 
entail the scalar MM or the pseudoscalar current M^M: 

q m = c M MMF^ v F uv , Of = C¥MMF^F UV , 

(2-2) 

= C^M^MF^F^ , = M^MF^F^ . 

For a Dirac fermion, the corresponding operators are again simply obtained by replacing M 
with D. Additionally, if DM is a Dirac fermion, we can write down an operator involving 
the tensor current Da^D, where a^ v = | (7^7^ — 7 i/ 7 M )- It reads 

Of = C^Da^DF/Fup . (2.3) 

In contrast to the case of scalar DM, the Wilson coefficients entering the definitions (2.2) 
and (2.3) are proportional to M~ 3 . A summary of all the effective operators considered in 
this paper together with their main properties is given in Table 1. 
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Table 1. Effective operators constructed out of DM and field strength bilinears together with 
their main properties. See Section 3 for a discussion of the di-photon annihilation cross sections and 
Section 5 for details on the direct detection cross sections. The cross sections for complex scalar 
(Dirac fermion) DM are smaller than the ones for real scalar (Majorana fermion) DM by a factor of 4. 

3 Annihilation cross sections 



In order to calculate the di-photon annihilation cross section of DM, we have to first compute 



the squared matrix elements for the process XX ~ > 77 with x 
operators introduced in the previous section, we find 



R, C, D, M. For the set of 
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4, F = D , 
16, F = M, 

F = D, 



(3.1) 



16, F = M. 
(M*) V (a -4mf,) , 



where s 



4ml + m 2 w 2 



■ ! ,,l x u re\ + ^( u rel) denotes the center of mass energy. In the following we will 
assume v TC \ ~ 1.3 • 10~ 3 c. We discuss the uncertainties concerning v re \ in Appendix A. 

In terms of the amplitudes (3.1) the velocity-averaged annihilation cross section can be 
written as 

<^rel> 77 = ^d£T > (3-2) 



647rm 2 



where the label z indicates the type of operator insertion and an additional factor of 1/2 
has been taken into account reflecting annihilation into two indistinguishable particles. The 
explicit results for the velocity-averaged annihilation cross section for each operator can be 
found in Table 1. Note that Dirac and complex scalar DM are not their own anti-particles 
resulting in an additional factor of 1 /2 when calculating the 7-ray flux from the annihilation 
rate, assuming equal amounts of DM particles and anti-particles. 
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Although we are ultimately interested in deriving bounds on the allowed di-photon 
annihilation signal imposed by direct detection, it seems worthwhile to estimate the scale 
M* suppressing the higher-dimensional operators. This can be done by setting the DM 
annihilation cross sections into 77 equal to the benchmark value ((7 x v ie \)^ ~ 10~ 27 cm 3 s , 
and then solving for the Wilson coefficients Cf. Since the higher-dimensional operators we 
consider can only be generated from loops of new particles charged under U (l) em > one obtains 
M* of around 100 GeV and 10 GeV for s- and p-wave annihilation, respectively. Note that this 
effective suppression scale does not directly correspond to the scale of new physics, because 
couplings as well as the multiplicities of new states have to be taken into account. The 
above numbers nevertheless suggest that scenarios leading to an observable 7-ray signal are 
difficult to construct. In fact, a model that gives rise to (a x v rc \)~ n ~ 10 _27 cm 3 s _1 should 
either contain a high multiplicity of states running in the loops that generate the effective 
operators or should feature resonant (s-channel) production in which case the mass scale M* 
gets replaced by the width of the mediator that is exchanged in the annihilation process. 
We will not discuss possible ultraviolet (UV) completions. 

4 Loop-induced effective interactions in direct detection 

The dimension 6 and 7 operators Of introduced in Section 2 lead to interactions between DM 
particles and nuclei which can be probed in direct detection experiments. There are in essence 
two kinds of effects which contribute to the DM-nucleus scattering cross section. The first one 
results from loop-induced couplings of DM to individual quarks and gluons, while the second 
one corresponds to Rayleigh scattering which stems from coherent interactions of the two 
photons with the entire nucleus [23]. As we explain below, the different contributions involve 
several widely separated energy scales, namely the high-energy scale M*, characterising the 
onset of non-standard dynamics, the heavy-quark thresholds itiq, the light-quark masses m q , 
the inverse nuclear coherence length Qq, and finally the momentum exchange q, involved in 
the low-energy scattering. In order to separate short- from long-distance physics, it will be 
useful to again employ an EFT. 

In the following we will present the basic steps of this EFT calculation, quoting explicit 
results for the operator O m only. Further details on the computation of the individual 
contributions itself are relegated to Appendix B. The generalisation to the case of the other 
operators is straightforward and will be presented at the end of the section. For the further 
discussion we also make the assumption that the scale of new physics M* is above the top- 
quark threshold mj. 

4.1 Operator mixing and threshold corrections 

The first important observation is that O m mixes under QED into the operator 



where m q denotes the mass of the SM quark q, which can have any flavor q = u,d, s, c, b, t. 
This mixing is determined by the UV pole of the Feynman diagram shown on the left-hand 
side in Figure 1. To leading logarithmic accuracy, one finds the following relation between 



Of = Cfm q MMqq 



(4.1) 



the Wilson coefficient of the operator of O q and that of O 




(4.2) 
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Figure 1. One-loop Feynman graphs showing the contributions to the DM-nucleus cross section 
in the case of O m . Mixing diagram generating Of (left), matching contribution giving rise to Oq 
(middle) , and matrix element describing the low-energy two-photon scattering of DM on the nucleus 
(right). See text for further details. 



where e q is the electric charge of the quark q and m% < [x < M* . Notice that we have assumed 
that the Wilson coefficient of Of vanishes at M*. 

We now evolve the Wilson coefficient Cf from M* down to mt, where we integrate out 
the top quark. Removing the heavy quark as an active degree of freedom gives rise to a finite 
threshold correction to the Wilson coefficient of the operator 

Of = Cf MMG a ^ u G a ^ u , (4.3) 

where G a '^ v denotes the field strength tensor of QCD. The relevant leading-order (LO) di- 
agram is shown in the middle of Figure 1. The corresponding matching is captured by the 
simple replacement [28] 

m t MMttCf(mt) -> MMG a '^G a ^Cf (m t ) , (4.4) 
with Cq given at next-to-leading order (NLO) by 

C G I (m t ) = -°^(l + S t )C^(m t ), (4.5) 

where 5t = lla s (mt) / ('In) [29]. Although 5t is formally of higher order, we will include such 
finite two-loop contributions in our analysis, because they are numerically non- negligible. 
Notice that once the top quark has been removed, the Wilson coefficient Cf and the corre- 
sponding logarithm is frozen at the threshold m< in the EFT. 

After the top quark has been integrated out, we then have to consider the mixing of 
the set of three operators O m , Of and Of. Like O m the operator Of mixes into Of. 
The relevant diagram is the QCD counterpart of the one displayed on the left in Figure 1 
with the photons replaced by gluons. As shown in Appendix B, the associated corrections 
are subleading and we will neglect them in what follows. The operator Of itself evolves like 
the QCD coupling constant, so that for scales mi, < [i < nit its Wilson coefficient takes the 
form 

CgM * ° ?M f (1 + 5 t ) In ( A 4) C»m . (4.6) 

At the scales mb and m c , the bottom and charm quarks are integrated out, which in 
full analogy to (4.5) results in finite matching corrections to Cf. Including all heavy-quark 
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threshold effects as well as the renormalisation group (RG) evolution of the individual Wilson 
coefficients, we find for 11 < m c the following expressions 

C M (fi) ~ C M (M*) , 



Cf^^-Se^lnf^jC^M). 



t V M / (4.7) 

Q V <9/ 

with g any of the light quarks u, d, s, while the sum over Q includes the heavy quarks c, b, t. 

4.2 Nuclear matrix elements 

Having determined the values of the Wilson coefficients at the hadronic scale [x = 0(1 GeV), 
we still have to evaluate the matrix elements of the operators O m , O^ 1 and Oq between 
nucleus states A. We begin our discussion with the matrix element (^4|O m |j4). The important 
observation in this context [23] is that in general the two virtual photons do not scatter off 
a single nucleon, but interact with the entire nucleus (see also [30] for a related discussion). 
This is illustrated by the Feynman diagram on the right in Figure 1. The strength of this 
interaction is to first order a measure of the total electric charge Ze of the nucleus, while 
effects related to its substructure (such as spin or magnetic moment) are suppressed by the 
nucleus mass. The matrix element of O m can hence be computed by treating the nucleus 
as a slowly moving charge source, applying well-known techniques of heavy-quark effective 
theory. The basic steps of the actual calculation can be found in Appendix A of [23] and will 
not be repeated here. To leading power in the nucleus mass and zeroth order in the velocity 
expansion, we find 1 

(A\O m \A) ~ 2^-aZ 2 Q Q F K&y (q)C M ^) , (4.8) 

where we have introduced the form factor 2 

l~2 f 1 f°° I 2 

F^ y (q) = -J- dx dl — — — 2 exp [-2 (l 2 - q 2 ((1 - x)x - §))] 

V TT Jo Jo [l + q (I — x)x) 

(4.9) 

I 2 — Q 2 (\ — X ) X _|_ I 

cosh (2lq (1 - 2s)) 2 - sinh (2lq (1 - 2x)) 

Iq (1 — 2x) 

normalised such that i ? R ay (0) = 1. Here q = \q\/Qo = q/Qo denotes the ratio of the three- 
momentum transfer q and the nuclear coherence scale Qq ~ 0.48(0.3 + 0.89^4 1//3 ) _1 GeV. 
For xenon with A ~ 131.3 one finds numerically Qq ~ 0.1 GeV. Notice that the coherent 
contribution (4.8) scales like Z 2 and is therefore more important for heavy nuclei. 

We now turn to the calculation of the matrix elements (A\Og \A). Following the stan- 
dard procedure in the DM literature, we evaluate these contributions at tree level. Neglecting 
small differences between the nuclear matrix elements of protons and neutrons, we obtain 

{A\Of\A) ~ 2m N A$ q F ne]m (q)CM ( M ) , (4.10) 



1 In (4.8), (4.10), and (4.12) the spinors associated to the DM bilinears are suppressed. 

2 Our results (4.8) and (4.9) differ slightly from the ones reported in [23]. The difference can be traced back 
to their equation (A-ll), where a factor of 2 is missing in the exponent after substituting (A-9) into (A-7). 
We also note that the sine function in (A-17) should be a hyperbolic sine as in our formula (4.9). 
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where — 0.939 GeV is the average nucleoli mass and A is the mass number of the nucleus, 
and the scalar form factor fJp is defined via mjqf^ q = (N\m q qq\N). Fn c \ m is the Helm form 
factor 

j? f — \ 3ji(qQ R,A) / i_2^,2.2\ (a 1 1 \ 

Faehniq) = — =77^ exp(-^q Q t ) , (4.11) 

qwo-tiA 

which depends on the effective nuclear radius Ra — \J (6.1 GeV" 1 A 1 / 3 ) 2 — 5t 2 and the nu- 
clear skin thickness t ~ 4.6 GeV -1 , while j\ denotes the spherical Bessel function of the first 
kind. Notice that a possible momentum dependence of the scalar form factors f^ q is not 
included in (4.10). Such effects together with other corrections arising at NLO have recently 
been computed in [31] employing the formalism of chiral perturbation theory (ChPT). While 
these corrections can be significant if there are large cancellations between the contributions 
from protons and neutrons, in our case their inclusion would have a minor impact only. The 
approximate result in equation (4.10) is hence sufficient for our purposes. 

The last missing ingredient in the calculation of the DM-nucleus cross section is the 
nuclear matrix element of the operator Oq. In contrast to (4.10), the first chiral corrections 
to (A\Oq\A) arise at the NNNLO level [31]. It is hence an excellent approximation to use 
the classic tree-level result of [28]. In our notation, one obtains 

(A\Og\A) * -2y^m N Af» G F Relm (q)C%(»), (4.12) 

where again the scalar form factors for protons and neutrons have been set equal. Following 
common practice, we have used the expression fJp G = 1 — fJp q for the gluon form factor. 

The nuclear matrix elements for the operator O r are identical to the results given 
above, while those for O c,D are smaller by a factor of 2 than (4.8), (4.10) and (4.12). The 
contributions from the operators O*, on the other hand, all vanish due to the anti-symmetry 
of the Levi-Civita tensor in F^. Finally, the pseudo-scalar operators O^J^ and the tensor 

operator Of generate only momentum-suppressed and spin-dependent interactions between 
quarks and DM at the one-loop level, respectively, which are irrelevant for direct detection. 

5 Direct detection cross sections 

We have seen in the previous section that only the dimension-6 operators O r,c and the 
dimension-7 operators O m,d lead to unsuppressed spin-independent (SI) interactions between 
DM particles and nuclei. For these cases the scattering cross section on the entire nucleus 
can be written as 

a *- n * n lf 1 \l, X = M ,D, (5 - 1} 

where \xa = ?T i J 4^x( m ^ + m x)~ l ls the reduced mass of the DM-nucleus system and n x = 2 
for x = -Rj M, while n x = 1 for x = C,D. The quantity f x takes the general form 

f x = L f{A\0*\A) + 5>4|0*|.A) + (A\O g \A)) , (5.2) 

q 

and includes the three contributions evaluated in the previous section. 
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Figure 2. Left panel: Comparison of the Rayleigh form factor (solid blue line) and the Helm 
form factor (dashed green line) as a function of q. The Rayleigh form factor shows a much stronger 
suppression towards large momentum transfer. The blue shaded region displays the typical values 
of q relevant for the XENON100 detector. Right panel: Recoil energy averaged form factors (F 2 ) 
or XENON100 as a function of the DM mass. The blue (green) band represents the result for the 
Rayleigh (Helm) form factor. See text for further details. 



Combining the relations (4.7) to (4.12) the following explicit expression for the effective 
coupling f x is found 



P 



o 



7T 



^Z 2 Q F Ray (q) - 3m N A (j + ^ (f» d + L N F Uelm (q) 



~j m N A$ G (-(A c L c + A t L t ) ■ ^ A/, Li, ) F u ,.uJ<]) 



(5.3) 



C X (M*) , 



where Lp = In (M 2 /mp) with P = N,c,b,t and Aq = 1 + Sq. Following [23], we will use 
a ~ 1/137 in our numerical analysis. The light-quark matrix elements can either be 
determined phenomenologically from baryon masses and meson-baryon scattering data or 
computed within lattice QCD (see [32] for a concise review). We adopt the values 



0.021, 0.038, 0.013. (5.4) 

The numbers for f^ u and f^ d have been obtained by averaging the values for the proton 
and neutron scalar form factors as given in [32], whereas has been taken from the recent 
lattice QCD study [33]. Compared to older values of typically f^ s — 0.14 [32] the number 
given in (5.4) is notably smaller. As discussed below, the resulting uncertainty associated to 
this choice is small. Notice also that in (5.3) we have identified the renormalisation scale // 
with the scale at which the matrix elements of the operators of Section 2 are matched onto 
ChPT or another low-energy effective Lagrangian describing the DM-nucleon interactions. 
For definiteness we have used the nucleon mass in this matching. In our EFT calculation 
we employ the MS scheme and correspondingly we take m c ~ 1.3 GeV, m& — 4.2 GeV and 
wit — 165 GeV for the heavy-quark masses. The associated threshold effects are A c ~ 1.35, 
A fe ~ 1.20 and A t ~ 1.10 for a s {m c ) ~ 0.399, a s (m b ) ~ 0.226 and a s (m t ) ~ 0.109. 

As an example, we calculate f x for interactions between the DM particles and xenon 
atoms. Using the parameters specified above, we find 

M 2 

0.57 + 0.09 In ( | r Relm{ 



e N 



f*~\ 1.69F Ray (g) 



(200 GeV) 2 



(q) } C x (M*)GeV. (5.5) 
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Figure 3. Left panel: Expected event rate in XENON100 assuming C x = l(T 6 GcV~ 3 . The blue 
(green) line indicated the Rayleigh (mixing) contribution, while the red line illustrates the actual 
event rate after interference of the two contributions. Contrary to most DM models, the event rate 
does not decrease towards larger recoil energies, giving a distinctive signature of this model. Right 
panel: Differential event rate for CDMS-II, measured in dru = events/ (kg day keV). The displayed 
results again correspond to the choice C x = 10~ 6 GeV~ 3 . 



We add that employing fx s — 0.14 instead of the value quoted in (5.4) would change the 
numerical coefficient 0.57 into 0.64, implying that the theoretical uncertainty in f x related 
to the strange-quark scalar form factor is about 10%. The errors on /iY and f^ d have even 
less impact on the obtained results. 

Our result for f x has some features worth noting. We first observe from (5.3) and (5.5) 
that there is a relative sign between the coherent contribution induced directly by O x and 
the corrections from Oq and Oq arising due to operator mixing and threshold corrections. 
The two types of effects hence interfere destructively. In fact, the extent of the interference 
depends sensitively on the value of the coherent form factor i ? R ay (g) compared to the standard 
form factor Fn e \ m (q). The g-dependence of both form factors is shown in Figure 2. From the 
curves it is evident that the coherent form factor is much more strongly suppressed for finite 
three-momentum transfer than the Helm form factor. 

Clearly, a realistic direct detection experiment will not be able to probe the scattering 
cross section at zero-momentum transfer. Instead, it will only be sensitive to nuclear recoils 
with energy = q 2 /{2mA) larger than the energy threshold £th of the detector. In order 
to make a realistic estimate of the relative importance of the two contributions we need to 
compare the two form factors at finite momentum transfer. In other words, we want to 
calculate the average value of the form factors, (F 2 ), in a given direct detection experiment. 
This goal can be achieved by calculating the recoil energy average of the form factors for 
different values of m x as outlined in Appendix C. 

The averaged Rayleigh and Helm form factors for XENON100 [34] as a function of 
m x are displayed in the right panel of Figure 2. We see that the averaged Rayleigh form 
factor is significantly smaller than its Helm form factor counterpart. The ratio of the two 
decreases from (F^ y ) / (F 2 elm ) ~ 0.3 for m x < m A to (F 2 ay )/(F 2 clm ) « 0.1 for m x > m A . 
To estimate the uncertainties of the averaged form factors, we have studied their dependence 
on the relative scintillation efficiency £ e ff in xenon and on the DM velocity distribution f(v). 
We find that the error related to C e g is the dominant individual source of uncertainty in {F 2 } 
and amounts to roughly 10% for large values of m x . The associated la (2a) error bands are 



- 10 - 




m x [GeV] m x [GeV] 

Figure 4. Left panel: Limits on C x as a function of m x from XENON100 (red) and CDMS-II 
(blue) for \ = M (solid) and x = D (dotted). The same limits apply to C x /m x for x = R (solid) and 
X = C (dotted). Right panel: Bound on (eri> rc i) 77 as a function of m x from XENON100 (red solid 
curve) compared to the current bound from Fermi LAT (black dashed curve) for Majorana DM. The 
projected sensitivity of XENON1T is indicated by the purple dotted curve. The green box shows the 
parameters for the claimed 7-ray signal from [4]. For Dirac DM, the bounds from XENON100 and 
XENON1T would be stronger by a factor of 2. All shown curves have been obtained by setting the 
new-physics scale M t entering the effective coupling f x equal to 200 GeV. 



displayed in light (dark) colours in the right panel of Figure 2. 

Because of the strong suppression of the Rayleigh form factor for finite momentum 
transfer, the two terms in (5.5) give comparable contributions to the differential event rate 
in a realistic detector. Typically, the contribution proportional to the Rayleigh form factor 
is larger close to the threshold, while the contribution proportional to the Helm form factor 
dominates at large recoil energies. Consequently, there will be large interference effects 
leading to a distinct recoil spectrum. This striking feature is illustrated in Figure 3 for 
both XENON100 [34] (left panel) and CDMS-II [35] (right panel). Note that the effect of 
the interference is much more pronounced in CDMS-II because of the much better energy 
resolution. Contrary to most DM models, we expect almost no events near the energy 
threshold, whereas the bulk of the signal is expected at relatively large momentum transfer. 

In spite of the large interference, we still obtain relevant constraints on the Wilson co- 
efficients C x both from XENON100 and CDMS-II using the maximum gap method [36]. The 
resulting bounds are shown in the left panel of Figure 4 and range from around 10 -3 GeV -3 
to almost 10 -6 GeV~ 3 . These constraints can be translated into bounds on the 77 annihila- 
tion cross section. In order to get a feeling for the resulting limits, we define <Tq = <7gj((7 2 = 0) 
and /o = f x (q 2 = 0)/C x . Combining (3.2) with (5.1), we then obtain 

4 5 4 F 

(0- Wrel) 77 < n S — , (<T V Ie l) J7 < ^--3-72 V TeU ( 5 - 6 ) 

jo z Jo 

where S = R,C and F = M, D. 

From the above relations it is readily seen that for scalar DM, corresponding to s- 
wave annihilation, the current sensitivity of direct detection searches is insufficient (by six to 
eight orders of magnitude) to compete with indirect detection experiments. In contrast, for 
fermionic DM which leads to p-wave annihilation, the present constraints from XENON100 
are, as a result of the factor in (5.6), comparable to the Fermi LAT constraints if 
m x < 40 GeV. For light DM, direct detection experiment hence start to indirectly probe 
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annihilation cross sections of almost (cx?; re i) 77 ~ 10 -29 cm 3 s _1 . This feature is illustrated by 
the red solid and the black dashed line in the right panel of Figure 4. Notice that towards 
larger masses, the bound becomes weaker in spite of the stronger limit on C x , because of the 
factor of m x appearing in (5.6). 

Nevertheless, if DM is fermionic and annihilates into 7-rays via the operator O x , the 
tentative 7-ray signal at m x = 130 GeV will be easily within reach of XENON1T [37] after a 
run-time of one year. This is indicated by the purple dotted line in the right panel of Figure 4. 
For this estimate, we have assumed that the XENON1T detector has identical properties as 
XENON100, except for the fiducial mass, which we take to be 1.1 tons. Assuming different 
detector properties close to the threshold does not affect the bounds for m x > 100 GeV. 
On the other hand, since we expect a significant number of events at large recoil energies, 
raising the upper bound of the DM search window to larger values of the primary scintillation 
signal SI would additionally increase the sensitivity of XENON1T. We have checked that 
the uncertainties in v re \, fjf and £ e g do not compromise our results. 

6 Discussions 

Motivated by recent results from searches for 7-ray lines and tentative claims of a signal, 
we have studied the direct detection cross section arising from interactions of DM with two 
photons. The di-photon annihilation cross section can be cast in terms of effective operators 
which allows for a translation into direct detection event rates. To this end we identified all 
effective operators giving rise to DM annihilation into two photons up to dimension 6 for real 
and complex scalar DM and up to dimension 7 for Dirac and Majorana DM and calculated the 
loop-induced scattering cross section in a model-independent way. Our operator analysis is 
valid under the assumption that the scale of new physics is sufficiently high and in particular 
that DM does not annihilate via an s-channel resonance. 

The loop-induced direct detection cross section receives contributions from two different 
types of effects. The first accounts for the fact that loop diagrams involving high-virtuality 
photons induce couplings of DM to quarks and gluons, while the second one corresponds to 
Rayleigh scattering that describes the low-energy interactions of the photons with the total 
electric charge of the nucleus. We discussed in detail the QED and QCD mixing of operators 
that leads to the former corrections, which are typically neglected in the direct detection 
literature. While the resulting contributions have the standard A 2 scaling, the Rayleigh 
contribution scales as Z 4 . 

At zero-momentum transfer the nuclear matrix elements for direct detection on xenon 
are comparable in size but have opposite signs, which leads to destructive interference. For 
finite momentum transfer, however, the form factors of the two contributions behave rather 
differently - the Rayleigh form factor gets suppressed much faster such that the overall event 
rates in xenon are dominated by the contribution due to mixing. Furthermore, the resulting 
recoil spectrum has a dip just at the lower end of the search window - a striking feature 
that, if observed, could confirm the interference of the two different form factors. For lighter 
targets Rayleigh scattering is more strongly suppressed due to the Z 4 dependence, while on 
the other hand the nuclear coherence scale Qq is larger which eases the form factor suppression 
for finite momentum transfer. 

The overall scale of the DM-nucleus cross section is proportional to the annihilation cross 
section into two photons and also depends on the leading partial wave of the annihilation 
process. For (<r x ?; re i) 77 ~ 10 _27 cm 3 s _1 , the induced direct detection cross section for s- 
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wave annihilation is too small to be observed in upcoming direct searches, even for spin- 
independent cross sections. Operators which lead to spin- or momentum-dependent cross 
sections are even less constrained by direct detection experiments. For p-wave annihilation, 
on the other hand, the DM nucleus cross section is just below the sensitivity of XENON100, 
but within reach of XENON1T. 
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A Velocity dispersion of Galactic DM 

Given a mass model of the Galaxy, consisting of a DM distribution p x (r) as well as baryonic 
contributions, the dispersion a xr = (vZ r ) of the DM radial velocity v X)T follows from the 
Jeans equation. In particular, assuming spherical symmetry and isotropy of the distribution 
function [38], 

<r(r) = ^l dr Px (r)-^>, (A.l) 

where $(r) is the gravitational potential generated by all mass components. The speed 
v x = yv^j, v X (j,^~ v x e * s distributed as a Maxwell-Boltzmann distribution proportional to 
v x exp [— v x /vq] with (v x ) = 3vq/2 and vq = 2a x>r (r). The dispersion of the relative velocity 
^rei = WxA ~ %.2| is simply (t> r 2 cl ) = ((w x ,i) 2 ) + {(v x ,2) 2 } ~ 2(^,1 • v x , 2 ) = 2(v x ) = 3ug, where 
we have used that (v Xt \ ■ v X) 2) must vanish by symmetry. For our analysis we adopt the 
benchmark value of vq = 220 km s -1 , which gives v re \ = y3vo ~ 1.3 • 10~ 3 c. 

In general, however, the velocity distribution of DM is a function of the galacto-centric 
radius r. For s-wave annihilation, one has (cr x v re i) = const, such that the r-dependence does 
not affect our results. For p-wave annihilation, however, one needs to estimate by how much 
the r-dependence of the velocity distribution enhances or suppresses the velocity-averaged 
annihilation cross section. 

To this end, we employed two state-of-the-art models [39] (model 1 and model 5) for 
the baryonic mass distribution of the Galaxy which have been shown to fit microlensing as 
well as rotation curve measurements. Model 1 consists of a bulge, a thin and a thick disk 
while model 5 consists of a bulge, a bar, a stellar disk as well as a gas disk. For both baryonic 
models we considered a Navarro-Frenk- White (NFW) profile [40] p x oc x~ a (l+x)~ 3+a (where 
x = r/r s , r s = 20kpc and a = 1) as well as an Einasto profile [41] p x oc exp [— 2{x a — l)/a] 
(where r s as above and a = 0.17) with the local DM density fixed to pr>M = 0.4GeVcm~ 3 . 
We compared (a x v re i)p x for the case of the full r-dependence to the case where only p x 
depends on r and (cr x v m i) is fixed to its local value. In the inner kiloparsec, the velocity 
dependence leads to a suppression proportional to r °- 5 ---°- 85 ) but when integrating over lines- 
of-sight and averaging over the target regions of [4] the contribution from the inner kiloparsecs 
is strongly suppressed. We find an overall enhancement of the velocity-averaged annihilation 
cross section with respect to the case with a fixed velocity distribution for an NFW (Einasto) 
profile of 15% (5%) for model 1 and 27% (30%) for model 5. This enhancement can be traced 
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back to the pronounced bulge/bar, especially in model 5, but barring the uncertainties of the 
mass modelling in the inner few kiloparsecs of the Galaxy, we chose to ignore the uncertainty 
related to this effect in our final results. 



B Mixing and matching of loop-induced DM operators 

In this appendix, we review in detail the mixing and matching of the following effective DM 
operators, 

QX = C* XX F^F^ , O* = qm qX XW , O x G = C^ X xG a ^G% . (B.l) 

Since by assumption the DM particle \ is a singlet under the SM gauge group, the \X parts 
of the above composite operators can be ignored for the further discussion. This feature 
implies that the RG evolution and the matching corrections to the Wilson coefficients in (B.l) 
resemble those well-known from Higgs physics. 

In fact, the renormalisation scale dependence of the Wilson coefficients is fully de- 
termined in terms of the QED beta function f3 e and the corresponding mass anomalous 
dimensions 7^ of the quarks 



n" + J 



— a^) = f3 e ^) = -2a^)Y J Pi n)fn{l ' ) 



n=0 



(B.2) 

t 00 / / NX n+1 



InmO) = 7 ^(/x) = - 5^7,,, 



^ n=0 v 

as well as their QCD counterparts (3 S and 7^. Note that for our definition of the beta 
functions, /3 e (/i) > while j3 s (p) < 0. Assuming that the Wilson coefficient Cq vanishes at 
the new-physics scale M*, one obtains for scales m t < fj, < M* the expressions 

w a(//)&(M,) 1 ; ' 
W = fa(M m ) - t£(m)) CHM*) + C*(M*) . 

These all-order results represent the QED analogs of the formulas given in [42] for the case 
of QCD. Here a, /3 e , and 7^ all correspond to quantities with the number of flavours set to 
those active at \x. Expanding (B.3) in powers of the electromagnetic coupling a and keeping 
only terms relevant for our analysis, one finds 

C*(^^U*(M«), 

9 ^ / ( \ \ (B - 4) 

C x( u) ~ 2 J^ 1 _ c x (M ) 

for the case C*(M*) = 0. Notice that the appearance of the rescaling factor a(/i)/a(M*) in 
the first expression guarantees the proper normalisation of the matrix element of O x at the 
scale [i. Using now the LO relations 

^ + (B.5) 
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and 7m°^ = 6ejj the results in (B.4) can be further simplified. In particular, neglecting small 
effects associated to the running of a, one arrives at (4.2). 

At each heavy-quark threshold tuq the Wilson coefficient C G receives finite matching 
corrections that are known to NNNLO [43, 44]. In the case of mt, one has including NNLO 
contributions [45, 46] 



a strut) 
12vr 



1 + 



lla s {m t ) , a 2 s {m t ) (2777 



4tt 



+ 



288?r 2 



m t 



(B-6) 



The same expression applies to the bottom- and charm-quark thresholds mj and m c after 
adjusting the number of active flavours rif. Since the second, flavour-dependent term in 
brackets represents only a relative correction of below 1%, we neglect such terms in our 
numerical analysis. Likewise, we do not include small finite threshold corrections to the 
Wilson coefficients Cq that start at order a 2 (see for example [42]). 

Once the top quark has been integrated out, one has to consider the mixing of the full 
set of operators (B.l). The solutions to the RG equations read in this case 

a(m t )/3 e (/i) 



C*(/i) 



a(n) p e (mt 



■C*(mt) 



ia s (m t 



a s (mt) j3 s {n) 
a s (n)l3 s (mt) 



{lm(m) -7m(/^))C X (m t ) 



(M))Cg(mt) 



(B.7) 



+ C}(m t ) 



where nib < fi < rUf. Performing again a series expansion in the coupling constants, using 
the QCD version of (B.5) and 7m°^ = 8, one derives from (B.7) the approximations 



C X (M*) 



-3e; 



a(M* 



In 



7T 



a(M«) a s (/x) 



7T 



7T 



1 + 



«(M,) / a B (m f ) Y 



7T 



In 



lla s (m t ) 



47T 



In 



7T 



77? 7 



In 



Mi 



C*(M*), (B.8) 



C X (M,) 



The appearance of the overall factor a s {n) in Cq is again crucial for the correct normalisation 
of the matrix element of Oq. Notice that the second term in the square bracket entering the 
expression for Cq is subleading and can be neglected for all practical purposes. Following the 
above discussion the implementation of the bottom- and charm-quark thresholds is straight- 
forward. Neglecting subleading 0(a 2 ) terms in (B.8) as well as the evolution of a, which is 
a very good approximation, the final expressions for the low-energy Wilson coefficients given 
in (4.7) are then readily obtained. 



C Direct detection for XENON 

The typical recoil energy in a given direct detection experiment is determined both by the 
expected recoil spectrum for DM interactions with the detector material and the properties 
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Figure 5. Left panel: Response function Res (E-r) for XENON100. The blue shaded region 
displays the acceptance window in the absence of fluctuations. Right panel: Normalised probability 
distribution of recoil energies p(Eji) for different values of the DM mass. 



of the detector itself. The recoil spectrum is proportional to the velocity integral 

9 m= r^^±^i, (CD 



where f(v) is the local DM velocity distribution evaluated in the Galactic rest frame, v& is the 
velocity of the Earth relative to this frame, and v = \v\. The minimum velocity required for 



a DM particle to transfer an energy Er to a nucleus is given by v m - m (E R ) = y tuaEr/ (2/i^). 

The performance of the detector, on the other hand, can be characterised by the response 
function Res (-Br.)- This function encodes the probability that a scattering process with 
recoil energy E R will lead to an observable signal within the DM search window of the 
experiment. For the XENON100 experiment, the search window is defined in terms of the 
primary scintillation signal SI by the condition that it should consist of between 3 and 20 
photoelectrons (pe) [34]. 

For a given recoil energy, the expected number of 5*1 photoelectrons is 

Sl(E R ) = E R L y C cS ^, (C.2) 

'-'ee 

where L y = 2.28pekeV~ 1 is the light yield of the detector, C e Q is the relative scintillation 
efficiency of liquid xenon, while 5 nr = 0.95 and S ee = 0.58 are the field-dependent quenching 
factors for nuclear and electron recoils, respectively. These values are taken from [34], while 
the best-fit curve for C e $ is taken from Figure 1 of [47] . Due to fluctuations, the actually ob- 
served number of photoelectrons in a collision may deviate significantly from the expectation 
value [48]. Assuming that the distribution of photoelectrons can be described with Poisson 
statistics, the probability to observe n photoelectrons in a collision with recoil energy E R is 
hence given by 

/foju--g<- a <*y> gl W'. (c.3) 

n! 

Consequently, the probability of a nuclear recoil with energy E R to give a primary scintillation 
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signal 51 within the search window is 

20 



Res(£ R ) = ^/(n,£ R ). (C.4) 

n=3 

The response function for XENON100 is shown in the left panel of Figure 5. For comparison 
also the acceptance window in the absence of fluctuations is indicated. 

It follows, that the observed distribution of recoil energies in the XENON100 detector 
is then given by the product 

p(E R ) = g(v min (E R )) Res (E R ) , (C.5) 

where the energy dependence of the cut acceptance in XENON100 has been neglected for 
simplicity. For light DM, g(v m i n (E R )) is a steeply falling function of E R , whereas the response 
function Res (E R ) is strongly suppressed for E R < E t ^ with the threshold energy determined 
by the requirement that S'l(E'th) = 3 pe. Consequently, the function p(E R ) will be strongly 
peaked at E R ~ E^ and as a result the typical momentum transfer for light DM amounts to 
q ~ ^2m A E th . 

For heavy DM, the function g(v m i n (E R )) decreases much more slowly and one needs 
a slightly more elaborate treatment. In order to deal with this case, we first define the 
normalised distribution of recoil energies 

p(E R ) = J^ R) . (C.6) 
Jo dE R p(E R ) 

Since v m i n (E R ) depends indirectly on the mass m x of the DM particle, so does p(E R ). The 
function p{E R ) is shown in the right panel of Figure 5 for different values of m x . One observes 
that for increasing DM mass the peak of the distribution is not only shifted to higher recoil 
energies, but that also the tail of the spectrum becomes much more pronounced. For an 
accurate description, one hence has to use momentum-averaged form factors 

(F)= / dE R p(E R )F(E R ). (C.7) 
J o 
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